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Abstract: We study various transitions in dilaton holography, including those associated 
with the spontaneous breaking of a global (superfluid case) or local (superconductor case) 
U(l) symmetry in diverse dimensions d. By analyzing the thermodynamics of the dilaton- 
gravity system we find that scale invariance is broken at low temperatures, as shown by a 
nontrivial hyperscaling violation exponent in the infrared; increasing the temperature we 
recover scale symmetry in a d dependent way: while for d = 2 + 1 a phase transition is found, 
for d = 3+1 the transition is rather a crossover. This is the expected behavior of QCD where 
the number of colors Nc equals three (although in our holographic calculations Nc — )• oo). 
When the U(l) is preserved and at low temperatures, the system is insulating for arbitrary 
d if the dilaton is appropriately coupled to the gauge field; for other couplings we also find 
a linear in temperature resistivity. We then determine the prediction of these models for 
several quantities in the superconducting phase: the DC and AC conductivity, the gap for 
charged excitations, the superfluid density, the vortex profiles, the coherence length, the 
penetration depth and the critical magnetic fields. We show that at low temperatures some 
of these quantities differ qualitatively compared with the corresponding models without the 
dilaton, although the superconductor is robustly of Type II. The ratio of the gap over the 
critical temperature of the superconductor is studied in detail varying d and the couplings 
of the dilaton and then compared with the BCS value. A holographic renormalization is 
required in d > 2 + 1 to compute some quantities (such as the AC conductivity and the 
penetration depth) and we explain in detail how to perform it. 
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1 Introduction and summary of the article 

The gauge/gravity duahty and its simplest reahzation, the anti de Sitter /conformal field 
theory (AdS/CFT) correspondence, have been extended in the last years to describe sys- 
tems at finite temperature T and density (associated with a U(l) symmetry). Much of the 
motivation for this setup comes from quantum chromo dynamics (QCD) at finite tempera- 
ture and baryon density, but more recent works have considered applications to condensed 
matter as well (for reviews see [1-7]). The virtue of this approach, which we will take 
advantage of in this paper, is that it allows us to have information on a strongly coupled 
gauge theory (in the limit of a large number of gauge degrees of freedom) by doing classical 
calculations in a gravity theory. 

The minimal bulk field content is given by the d + 1 dimensional metric and gauge 
field, which are dual to the d dimensional energy momentum tensor and U(l) current. 
An interesting extension, however, is provided by adding a real bulk scalar 4> (the dilaton), 
which couples to the gauge field essentially by changing its gauge coupling: g — t- g{(p). Even 
without modifying the ultraviolet (UV) properties of the theory, the dilaton can modify the 
infrared (IR) limit relaxing the (often too strong) constraints that conformal invariance 
imposes: a nontrivial (thus nonconstant) dilaton corresponds indeed to a nonvanishing 
beta function in the dual field theory. Most notably the presence of (j) can allow a charged 
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asymptotically AdS black hole (BH) to have vanishing entropy in the low temperature limit 
[8-10], contrary to what happens with the minimal field content. Dilatons, moreover, are 
typically present in the spectrum of string theories, which are natural UV completions for 
such effective holographic models. 

In this paper we are interested in various transitions occurring in dilaton holographic 
models in generic space-time dimension, including those that are associated with the spon- 
taneous breaking of the U(l) symmetry. To this end we introduce a charged scalar dual 
to a condensing operator. If the U(l) current is identified with the electromagnetic one, 
such breaking corresponds to superconductivity. The main motivation beyond holographic 
superconductors [11-13] is that they might lead to a better understanding of high-T su- 
perconductors, which certainly evade the weak coupling assumption of the theory proposed 
by Bardeen, Cooper, and Schrieffer (BCS). Cuprate high-T superconductors are obtained 
by doping so called Mott insulators and in this process one typically observes an insula- 
tor/superconductor transition [14]. In ref. [15] we showed^ for d = 2 -\- 1 that an (ap- 
proximate) insulator /superconductor transition can be obtained in dilatonic holographic 
superconductors for small enough T and if g{(f)) is big enough in the IR. Here we generalize 
the results of [15] to various space-time dimensions, including the physically interesting case 
d = 3 -|- 1, and obtain the prediction of holography for important observables (such as the 
energy gap for charged excitations in the conductivity and the coherence length), which 
were not computed in [15]. 

Another necessary ingredient for a realistic description of superconductivity through 
holography is the presence of a dynamical electromagnetic (EM) field. Dynamical gauge 
fields can be included by changing the usual Dirichlet AdS-boundary condition for the 
corresponding bulk field (suitable for global U(l) symmetries) with one of the Neumann 
type. Ref. [18] explained how to introduce a dynamical EM field in the original holographic 
superconductor model [11, 12]. Here we make use of this method to define and study 
genuine superconductors in dilaton holography for various space-time dimensions. We also 
compare the resulting theory with the one with standard Dirichlet boundary conditions, 
which corresponds physically to a superfluid rather than to a superconductor. Superfluids 
indeed, in a sense that will be specified below, are equivalent to superconductors in the 
limit in which the EM field becomes nondynamicaP. 

Let us give in detail the results of the present paper (and provide at the same time the 
organization) . In the following section we define the holographic model and the gauge/gravity 
duality; we also recall the method of [18] to introduce a dynamical gauge field. As we will 
see, for d > 2 + 1 this will force us to introduce a kinetic term localized on the AdS-boundary, 
which absorbs UV divergences through a holographic renormalization; more generally, such 
term is needed whenever the EM field has a nontrivial dependence on the space-time point, 
even if the EM field is nondynamical: for example this is the case when we determine the 
conductivity by computing the response of the system under a fixed (thus nondynamical) 

^Other studies of holographic superconductors with dilatons in d = 2 + 1 can be found in [16] (see also 
[17] for a previous study, which included the broken phase only). 

^Therefore, in order not to duplicate symbols we will adopt the terminology of superconductivity rather 
than that of superfluidity, unless otherwise stated (see however table 1 for a dictionary between the two). 
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EM wave. In the same section we also study the solutions and the thermodynamics of the 
dilaton-gravity system when the gauge field and the charged scalar are set to zero. We 
consider the general asymptotically AdS static BHs with d — 1 rotation and translation 
symmetry, which were recently presented in [19, 20]. Fixed the parameters of the model 
there is only one solution with active dilaton (dilaton-BH henceforth), but the fact that 
the space is asymptotically AdS implies the usual AdS Schwarzschild black hole (S-BH) is 
always a solution. We therefore compare the free energies of the S-BH and dilaton-BH to 
determine the energetically favorable phase. We find that at low temperatures the system 
prefers the dilaton-BH and in so doing breaks conformal invariance as shown by the emer- 
gence of a nontrivial hyperscaling violation exponent in the IR. Increasing the temperature 
we observe a physically different behavior for d = 2 + 1 and d = 3 + 1: while in d = 2 -|- 1 
there is a phase transition and the system jumps to the S-BH at some critical tempera- 
ture, for d = 3 -|- 1 we rather have a smooth crossover from the dilaton- to the S-BH. It 
is intriguing to notice that the confined/deconfined transition in real world {Nc = 3) QCD 
is expected to be a crossover. However, our results are in the limit of a large number of 
gauge degrees of freedom, which suggests that our dual theories cannot be identified with 
the large Nc extrapolation of QCD. 

In section 4 we study the simplest example of holographic superfluid phase on top 
of the dilaton-BH: a homogeneous configuration with condensed charged operator at finite 
temperature and chemical potential fi. As we will see a finite U(l) density (// 7^ 0) is needed 
in order to have a regular solution of the field equations. This configuration always exists 
when < T < Tc, where Tc is the critical temperature of the superconductor. Moreover, in 
this window of temperatures it is energetically favorable with respect to the normal phase 
(where the vacuum expectation value of the charged operator vanishes) like for the S-BH 
[11, 12]. However, we observe that the dilaton has an impact on the condensate, which 
depends on its couplings to the other fields and on d. 

We then analyze the DC and AC conductivity on top of these BHs in section 5, consid- 
ering both the normal and broken phase. While the presence of the horizon implies that the 
normal phase is conducting, if the coupling of the dilaton to the gauge field is appropriately 
chosen we have an (approximate) insulating behavior for the dilaton-BH. This effect is even 
stronger for small T and is enhanced in d = 3 -|- 1 with respect to the d = 2 + 1 case. For 
other choices of the dilaton couplings one can also realize a linear in temperature resistivity 
for small T. In the broken phase we analyze the superfluid density Ug and study in detail 
the gap of charge carriers ojg. In particular, we compute the ratio ojg/Tc to compare it 
with the BCS value, 3.5. We find, like in the S-BH [11, 12], a prediction which is robustly 
bigger than 3.5, but the dilaton does have an impact on such ratio. In particular, when the 
couplings of the dilaton are such that the normal phase is (approximately) insulating ujg/Tc 
is enhanced compared to the value found in the S-BH; there is also a range of parameters 
for which it is instead suppressed. As we have already pointed out, the calculation of the 
frequency dependent conductivity in d > 2 -|- 1 requires a holographic renormalization. We 
describe in detail how to perform it and show that the DC conductivity, Ug and ujg are 
renormalization scheme independent. 

Finally, in section 6 we study the response of the dilaton holographic superconductors 



-3- 



under an external magnetic field H and compare it with the case without the dilaton 
[13, 18, 21-28]. On the one hand, the inclusion of magnetic fields is important as some 
of the most famous properties of superconductors emerge when H ^ 0. On the other 
hand, the inhomogeneous configuration which are triggered by H allow us to illustrate the 
main differences between the global (superfluid) and local (superconductor) case. As we 
mentioned, in order to do this for d > 2 + 1 we have to properly renormalize the EM 
field at least in the superconductor case. We find the Meissner effect and study the vortex 
phase including nontrivial couplings of the dilaton to the gauge field and charged scalar. We 
obtain the coherence length penetration depth A' and the critical magnetic fields Hd and 
Hc2- The latter observables allows us to show that also these holographic superconductors 
are of Type II, both for d = 2 + 1 and d = 3 + 1, like those studied in [18, 27]. However, 
we find that A' and Hd have a qualitatively different T — )• limit with respect to the S-BH 
[27]: A' may remain finite if the couplings of the dilaton are appropriately chosen both in 
d = 2 + 1 and d = 3 + 1 and Hd shows an analogue difference. Interestingly, for these 
values of the couplings the normal phase is (approximately) insulating. 
Section 7 contains some outlook for future work. 



2 The holographic model 

We consider the standard Einstein-dilaton action coupled to a U(l) gauge field and a 
charged scalar ^' in d + 1 space-time dimensions: 

s = Jd^-^.V^[^[n-id.^f-vi^)] -^^l,-'0\D.M'], (2.1) 

where TZ is the Ricci scalar, Gn is the Newton constant, </) is a real scalar (the dilaton) 
and we introduced J^ap = daAp — djsAa and = da — iA^- A negative cosmological 
constant A is included in 1^(0), and the AdS radius is defined by A = —d{d — 1)/L'^. We 
define (p such that V{0) = A. Also the functions Za{4>) and Z^[(f)) are generic and we 
choose the normalization of the matter fields in a way that Z^(0) = .Z'^(O) = 1. For the 
Zs we only require to be regular and nonvanishing for any in order for the semiclassical 
approximation to be valid. For simplicity, we have not added any potential for ^. 

In this article we focus on the most general static asymptotically AdS planar black hole 
with d — 1 dimensional rotation and translation invariance, which has the following form: 

ds" = W{z) i-fiz)df + V + — \ , = ^(z) , (2.2) 

where dy'^ = dijdy'^dy^ and z is the holographic coordinate. The AdS-boundary is located 
at z = 0: we have W{z) ~ L'^/z'^, f{z) ~ 1 and (j){z) ~ for z ~ 0. 

2.1 The gauge/gravity correspondence 

Such d + 1 dimensional classical gravitational theory is supposed to be dual (through the 
gauge/gravity correspondence) to a d dimensional strongly coupled CFT with a large num- 
ber of gauge degrees of freedom. In the CFT side, the sources which generate the Green's 
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functions are identified with the AdS-boundary values of bulk fields; the generating func- 
tional itself is conjectured to be the gravity action (2.1) computed on a solution to the field 
equations with such AdS-boundary values as boundary conditions. For example, 

^\z=o = *o , A^\,=o = ' (2-3) 

can be identified with the sources of a charged operator O and the current respectively. 
One then obtains the following vacuum expectation values: 

W = ^^'-'-F.^U=o , {O) = ^z'-^D,^U=o . (2.4) 
5 9 

In this formulation of the gauge/gravity correspondence the boundary gauge field is 
an external source for the current and as such is a nondynamical field, introduced through 
a Dirichlet boundary condition (the second equation in (2.3)). However, it is possible to 
make dynamical by substituting this condition with one of the Neumann type. 



L 



d-3 



9' 



2=0 ef 



+ JLt = , (2.5) 



where is an external nondynamical current and we have added a bare boundary kinetic 



term ^8^^''^' 



for generality [18]. Condition (2.5) requires to add 



1 ~2 



4e 



(2.6) 



to (2.1). In the case d = 2 -|- 1 the bare boundary kinetic term is not necessary [29] and one 
can preserve conformal symmetry in the UV; in such case the dynamical is an emergent 
gauge field [18]. When d = 3 -|- 1 the situation is different. One way to see it is to notice 
that the current ( J^) contains a logarithmically divergent piece, 

^dzAfj^\z=o = -d''T,yf^\nz\z=o+ finite terms . (2.7) 
This divergence, however, can be reabsorbed in ef,: 

-IT = ^ -\ — ^ln2;L=o + finite terms , (2.8) 
Cfe 

where cq is the gauge coupling in the normal (^ = 0) phase. In this case behaves very 
much as the four dimensional photon we are used to: the logarithmic divergence in (2.8) 
corresponds in the dual field theory to the logarithmic divergence of the photon self-energy 
at very high momenta. For generic d > 3 -|- 1 the situation is similar to the d = 3 -|- 1 case, 
although the logarithmic divergence is substituted by a stronger one. 

This method to render dynamical has been originally applied to the holographic 
superconductor model [11] (including diverse dimensions [12]) in [18] and we refer to this 
paper for further details. Then, it has been applied to holographic p-wave [30-33] and 
gapped [34, 35] superconductors in [36, 37] and [27] respectively. In this paper we will use 
it to describe genuine superconductors in dilaton-gravity, extending the d = 2 -|- 1 analysis 
of [15] to various dimensions and to compute additional observables. 
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2.2 Solutions of the dilaton-gravity model 

We will work in the limit Gat — )• such that the effect of and ^ on Gap and (p can be 
neglected. In this limit the Einstein and dilaton equations are 



n. 



G. 



■al3 



[n - (a 



.^V^; ' - V{<P)\ - da<P dp<p = 



and 



-G 



-GG' 



dV 



(2.9) 



(2.10) 



respectively and the most general asymptotically AdS black hole of the form (2.2) solving 
these equations can be written as 



W{z) 



(1 + zjL) 

[(i+^/^r-i] 



l + 3(- 



2 ' 



+ (l + ^/L)^ 



1 



ln(l + zjV) , 

{l + z/LY {l + z/L)- 



for d = 3 [19] , while, for d = 4, we have [20] 



W{z) 



v\l + z/L) 
[{l + z/LY -l\ 

1 - 24 



3(z^2 _ 1) 



u + 2 



ln(l + z/L) 



,(2.11) 



16 



1/2 - 25)(9i/2 - 25) 



+ 



(1 + z/Lfl^ 



(1 + z/L)3'^/2 

3(3z/ + 5) 



+ 



(l + z/L)-3^/2 (l + z/L)'^/2 (l + z/L)-''/2' 

3(3i/ - 5) 1^ + 5 v-b 



(2.12) 



where > 1 is a dimensionless parameter and > has the dimension of a length. 

As usual we define the black hole horizon Zh hy f{zh) = In order to study the 
theory at finite T, we perform the Euclidean continuation to imaginary time it with period 
(3 = 1/T. Then T = |/'(z/j)|/47r, where the prime denotes the derivative with respect to z. 

In order for the functions W{z)^ f{z) and (piz) to be solutions to the field equations 
the potential has to be appropriately chosen: we have an exponential form^ 



(2.13) 



■t=i 



where V^'^\ 6^'^^ and the natural number Nd are constants. For d = 3 we have = 6 and 



L 



1-1/ 

> 

*2 + 1/ 



12- 



, (1 + ^^) 



l/2(i/-2) 



U+l ,2 3z/2_(j.2_4) 3/^3 



3Z/2 



4) /L3 



z/2(i. + 2) 



, (2.14) 



^Since the potential is even under (/> — >■ —0 we also have the solution in which we change the sign of 4> in 
(2.11) and (2.12). However, without loss of generality we consider only the sign choice of (2.11) and (2.12). 
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where the entries on the right hand side correspond to i = 1 , . . . , A'^^ respectively and 

-I/ - 1 



- 6P 



-1, 1 — V , u — 1 , 1, V + 1 . 



(2.15) 



For d = 4 we have instead = 7, 

{i)_da{u+l) da{v-l) f>da{v'^-l) 



db 



db 



2(3z^ - 5) ' 2(3z^ + 5) ' 9z^2 _ 25 

dc{y - 1) dc{v + l) 

(i/ + 5)(3i/ + 5) ' (jv-5)(3i/-5) ' 



3z^ + 5 



3z/ 



where 
da = 
and 



9z/2 - 25 



+ 



4(12LVz, 



2 /^2\2 



2 (z^2 _25)(9i/2 _25) 



(12L)2 5(1 - ;^2) 
1/3 4 zy2 _ 25 



(12L)2 



3(z.2 



1 - , V + 1 , 1 , 

-(3 + i/)/2, (i/-3)/2, -3(i/ + l)/2, 3(z^-l)/2. (2.16) 

We see that for i/ = 1 we recover the S-BH: W{z) = L? / z^, f{z) = 1 — (z/zg)^ and 
(j) = 0. Therefore u — 1 measures how strongly scale invariance is broken in the IR. Indeed, 
if we take the zero temperature limit we find f{z) — )• 1 and Ty(z) in the deep IR, z ^ L, 



can be approximated by Ty(z) 



^'{L/zY~^^, showing that v — 1 ^ corresponds to 



nontrivial hyperscaling violation exponent 9 [38-40]. The explicit expression of 9 is 

9 = ^{d-l) , 



(2.17) 



showing that 9 is d-dependent and always 9 < 0. Hyperscaling violation exponents with 
these properties can be realized in string theory [40]. 

3 Transitions in dilaton-gravity 

Since V{0) = A < 0, in addition to the solutions with nonvanishing dilaton in eqs. (2.11) 
and (2.12), we always have the usual S-BH, independently on the value of z/. In this section 
we determine which solution among the S-BH and the dilaton-BH is energetically favorable. 

In order to do so we have to study the free energy F of the system, which is given by 
F = TSe, where Se is the on-shell Euclidean action. By using (2.1) we have 



1 



d'^x. 



-n + {do,<j)f + v{(t)) 



(3.1) 



where the integral is performed over the spatial dimensions. By taking the trace of the 
Einstein equations in (2.9), we obtain an expression for 7^, which inserted in (3.1) gives 



By substituting the ansatz in (2.2) into this equation we then obtain 
^ Vd-i l"<^ 



(3.3) 



where Vd~i is the volume of the d — 1 dimensional space of the CFT and the parameter e 
represents a UV cut-off. Let us start with the d = 3 case for which the integrand in (3.3) 
is —W'^{z)V{4>{z)) and close to the AdS-boundary we have 



W'\z)V{^{z)) ~ ^ + ^^Y20l""'^^ + <^(^) = 3 dilaton-BH) 



t/-^+19-20i/^ 



-W'^{z)V{(t){z)) ~ 



6L2 



(ci = 3 S-BH) . 



From this equation we see not only that the UV asymptotics of the two solutions are the 
same but also that the UV divergences cancel in the difference F^ — Fs^ where F^ and Fg are 
the free energy densities for the dilaton-BH and the S-BH respectively. Performing explicitly 
the integral in (3.3), we find that the S-BH is energetically favorable at high temperatures 
while, decreasing the temperature below a critical threshold, a phase transition occurs such 
that the system jumps to the dilaton-BH. This situation is described in figure 1. 



L(Fj-F,) 0.00 



phase transition 

d=2 + l 



V = 1.2 



0.00 0.05 0.10 0.15 0.20 0.25 



Figure 1. The difference between the dilaton-BH and the S-BH free energy densities (up to the 
gravitational constant I/IGttGn) vs. the temperature for d — 2 + 1. 

When d = 4 we have a qualitatively different physics. The integrand in (3.3) is 
-2W^^^{z)V{4>{z))/3 and close to the AdS-boundary we have 

-IW'/Hz)V{m) + + ^'-"^'^illi''''^ + 0{z) {d = 4 dilaton-BH) , 

-lW^I^{z)V{(t){z)) = ^ (d = 4 S-BH) . 

Although the dilaton- and S-BH are both asymptotically AdS and as such their free energy 
densities have the same leading UV divergence, the difference F^ — Fg is not UV finite 
for d = 4, but is dominated by the next-to-leading divergence of the dilaton free energy, 
(1 — v^^je^. Since 1 — z^^ < for the dilaton-BH, it follows that the solution with active 
dilaton is always energetically favorable with respect to the S-BH. Nevertheless, when the 
temperature is increased, Zh becomes small in units of L and the dilaton-BH resembles 
more and more the S-BH. In other words we observe a crossover from a situation where 
scale invariance is broken to one where scale invariance is gradually recovered. The latter 
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region of the phase diagram can be interpreted as a deconfined phase. The fact that here we 
have a crossover rather than a phase transition is intriguing as deconfinement is expected 
to correspond to a crossover in the phase diagram of real world QCD (with the number of 
colors Nc = 3) [41]. However, we should point out that increasing Nc in QCD one observes 
a stronger deconfining phase transition, instead of a crossover [42-44]. At the very least 
this suggests that the dual of the holographic models we are studying is not the large Nc 
extrapolation of QCD, but of distinct field theories. 



4 Homogeneous superfluid transitions 

Let us now look for the simplest superconducting ^ 0) solutions and study how and 
when phase transitions from the normal phase to such superconducting solution occur. 

As we will see, there is no regular solution with 7^ if the gauge field is set to zero 
and a way to avoid this problem is to have Aq ^ 0. Therefore, we take the ansatz 

^ = ^P{z), Ao = Ao{z), (4.1) 

where is a real function, and the spatial components of A^ are set to zero. The require- 
ment that AaA'^ is finite implies ^oU=zh = (see ref. [11]). Thus, in order to keep Aq 
nonvanishing somewhere we impose ^0 7^ at the AdS-boundary; remembering also that a 
superconducting solution breaks the U(l) spontaneously rather than explicitly, we impose 

^-0 = , ao = Ai , (4.2) 

where ^ is a nonvanishing chemical potential. So the system is at finite density. 
The ansatz in (4.1) corresponds to the field equations 

a4VF('^-i)/2/Z^(0)9,V.) Al ^ _ {W^''-^)I^ZA{mA,) 2W^IJ^ _ 

and here we can see that there are no regular superconducting solutions with Aa = 0; 
indeed setting ^0 = in the equations above we find oc (VK^'^~"^^/^/.^^((^)) , which is 
singular at the black hole horizon. 

The boundary conditions in (4.2) are not sufficient to determine the solution to eqs. 
(4.3). However, the regularity at the black hole horizon provides us with other two condi- 
tions. Indeed, from the second equation in (4.3) we see that in order for Aq (and its first 
and second derivatives) to be regular at z = Zh we should impose = 0, as previously 
found by demanding AaA"' to be regular. More precisely we see that Aq ^ as z ^ Zh 
at least as fast as /. Then, we can see that for z ^ Zh the first equation in (4.3) reduces 
to W^'^^'^^/'^fZ^{4>)dz'il^ = with corrections which go to zero as z ^ Zh. Since W, f and 
Z^{(l)) are also regular at the black hole horizon, we find 

d,4^ = 0, Ao = {atz = Zh) . (4.4) 

It is worth noting that while (4.2) represent external conditions and as such are chosen 
according to the physics we want to describe, (4.4) are rather fixed conditions (as required 
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Figure 2. The condensate as a function of ^L. We set — 1.2, ZA{(t>) = e'^'^ and Z.ff, — e""^. We 
choose T ~ 0.010/L in the inset and T ~ 0.15/L everywhere else. 



by regularity). From the dual CFT perspective, a-iid are external sources while the 
conditions in (4.4) (and eqs (4.3)) correspond to the dynamical response of the CFT. 

We have solved eqs. (4.3) with the boundary conditions (4.2) and (4.4) and we plot the 
resulting (O) as a function of in figure 2. We see that, like in the holographic models 
without the dilaton [11, 12], there is a critical value of yuL above which the condensate 
becomes nonzero. Such critical value, however, is a function not only of T/L, but also of u, 
Za and and may be inverted to have Tc/ L in terms of the chemical potential, v and the 
dilaton couplings to and ^. Moreover, we have checked that, when the superconducting 
solution exists, it is energetically favorable with respect to the normal phase, as it can be 
shown by comparing their free energies. From figure 2 we see that the effect of the dilaton 
functions couplings is substantially different in the d = 2 + 1 and d = 3 + 1. 



5 Conductivity 

We now turn to the (complex) conductivity a. We will study this quantity both in the 
normal and superconducting phases. The conductivity in the normal phase will be given 
explicitly in section 5.1; this will allow us, among other things, to investigate the nature 
of the superfluid phase transition studied in the previous section. The conductivity in the 
broken phase will then be presented in section 5.2; from this quantity one can indeed observe 
an infinite DC conductivity and that the system develops a gap for charged carriers. 
Consider a small time-dependent gauge field along a spatial dimension x, 

A,{t,z)=A{z)e'^^P'-'^-'\ (5.1) 

which is induced by a small EM field at the AdS-boundary, ax{t). Here A and p are real 
functions of z. The system responds creating a current linear in Ux- {Jx) = o'Ex, where 
Ex = —dtax is the electric field, a can be computed once the equations of the other 
fields (other than the EM one) are solved on the EM background, ax{t); therefore, the 
conductivity does not depend on the dynamics of the EM field. 

The current (Jx) can be computed with the gauge/gravity correspondence (the first 
equation in (2.4)), to obtain 



L 



d-3 



A ^ ■ A '\ 

Aa; — + luAxP I 



(5.2) 
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having chosen a gauge such that = 0. Notice that this expression is generically UV 
divergent when d > 2 + 1 because the EM field we are considering is time-dependent: by 
looking at eq. (2.7), which holds for d = 3+1, we see that (J^) is logarithmically divergent 
whenever d^J-^^ 7^ 0; that singularity is replaced by an even stronger one when d > 3 + 1. 
We remove this divergence by adding a boundary counter term of the form given in the 
first term of (2.6). In the particular case d = 3 + 1 this means that l/e;, is logarithmically 
divergent in order for eo (defined in (2.8)) to be finite. This is nothing but a renormalization 
procedure and make expressions like (5.2) finite. In the case of {Jx) we have 



d-3 



^23^ + iujAxP 



A 



(5.3) 



By substituting eq. (5.3) in the definition of the conductivity a = {Jx)/Ex we obtain 



^Z-d 



9' 



P 



(5.4) 



z=0 



So the real and imaginary parts of the conductivity are given by 



Re [a] 



L 



d--i 



z=0 



lm[a] 



L 



d-3 



A' 



3-d -^.^ 



(5.5) 



z=0 



From this result we see that the boundary kinetic term only changes Im[cj] and it does 
so only by adding a term proportional to u. This implies that physical quantities such 
as the DC conductivity, the superfluid density Ug and the gap of charge carriers (in the 
superfluid phase) are independent on how we choose l/e^, namely they are independent on 
the particular renormalization scheme. 

Since ax{t) is small we can determine all the quantities needed to compute a in (5.4) 
by solving the linearized Maxwell equation 



dz {W^''-^)l^fZA{<l^)d,Ax] 



A 

+ j^Ax 



. 



(5.6) 



In the last term of the equation above we are taking into account the possibility that the 
system is superconducting, V 7^ 0. Eq. (5.6) has to be solved with appropriate boundary 
conditions. At the black hole horizon we impose, as required by regularity [1], that the EM 
wave is going towards the horizon (ingoing boundary condition), that is 



P{z)f{z) 



1 



as 



(5.7) 



Inserting (5.1) into (5.6) (and after that dividing by Ax) we obtain an equation whose 
imaginary part can be solved analytically, to obtain 



p'{z) = 5'ReH 



M(o)V 



W(<i-^)/^z)ZA{m)f{z) \A{z)J ' 



(5.8) 



where we used the expression for Re[cr] given in (5.5). A consequence is that Re[c7] cannot be 
zero without violating eq. (5.7). Therefore, the black hole horizon generates an imaginary 
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part of the amphtude A{z)e'^'^^^^\ which in turn results in a nonvanishing Re[(T]. In the 
fluid mechanical interpretation Ke[a] 7^ corresponds to a fluid, as opposed to a solid which 
has a vanishing Re[a] (see [27]). Inserting (5.8) into (5.6) we obtain an equation for A, 



where 

s(^)^ln4^-iln(p'M^(^-3)/2 



A{0) 2 V Jz=o' 
which, once inserted in (5.8), gives an expression of p' in terms of s, 

P'= W('i-^)/^ZA{(l^)f ■ ^^'^^^ 
We solve eq. (5.9) with the regularity boundary conditions (at z = z^) 

s- ^lnrVF(^-3)/2z,^ . _ 2fWZ^i^yL' - [2Z'^ + jd - 3)ZaW'/W] 

s---ln^W Za), s- Z^(4a;2 + f2) ■^^■^^> 

The first of these conditions comes from the requirement that the EM wave must go towards 
the horizon, eq. (5.7), and eq. (5.10); it ensures that the terms going like l/{z — z^)^ in eq. 
(5.9) vanish. The second one ensures that the next to leading terms, ~ l/(z — z/^), vanish 
as well. From the solutions to (5.9) and (5.11) we can compute a with 



1 / L'^~'^ e'"^^ \ 



1 /Z\3-rf , oj' 
b 



(5.12) 



2=0 

Before presenting a explicitly let us observe that the behavior of Re[(T] at high frequen 
cies^ can be easily determined. Indeed, in this limit (5.9) and (5.11) tell us 

1 , / ry2 



^hi[W''-'Z\'j+0{l/uj') (5.13) 

so one easily finds Re[(T] — >• 1/5^ when d = 2 + 1 and Re[o"] — t- +00 for c? > 2 + 1. This result 
is in agreement with [12], where they analyzed a in the superconductive case (and in the 
absence of the dilaton). These findings, however, do not rely on the fact that the system 
is superconducting and hold in the normal phase as well: as expected the system does not 
know if the symmetry is broken at high enough energies. 

5.1 Conductivity in the unbroken phase 

Let us now focus on the normal phase (see also the related study [45] for a charged BH). A 
particularly interesting case is the low frequency limit a; — )• as it allows us, in particular, 
to obtain the DC conductivity. We find 

1 



lim Re[o"l — „ 



^(d-3)/2^^ , limIm[cT]=0. (5.14) 



*Im[a"] at high frequencies is instead renormalization scheme dependent as different choices for the finite 
part of l/ei, can change the term proportional to uj. When d < 3 + 1 and the renormalization is not required, 
taking e;, — >■ 00 (that is without putting the boundary kinetic term) one obtains \m[a\ — > 0. 
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Figure 3. The conductivity as a function of the frequency in the unbroken phase. (O) — 0. We 
set T ~ 0.010/i, V = 1.2, = e''^ and = 1. Moreover, we choose jel = L ln(z/i)|^=o 

for d = 3 + 1, while Cb/g — > oo for rf = 2 + 1. 



The result for the DC conductivity in (5.14) is interesting as it tells us that this quantity 
is suppressed when Zj^ is small at the black hole horizon, resembling an insulator, while it 
is enhanced if Za\z=zk is large. 

These effects become strong in the low temperature region where the dilaton is log- 
arithmically large, eqs. (2.11) and (2.12). Indeed, taking the typical exponential form 
Za{(\>) = e'*"^, where 7 is a constant, we find at low temperatures 

lim Re[a] ~ T'^^JW^^ = 2+I) , lim Re^ ~ t^^+^^/^-t^/sF^^ = 3^1) _ 

Therefore, taking 7 < (that is suppressing Z^|2=2^), we have an almost insulating ma- 
terial. The d = 2 -|- 1 result reproduces that found in [15]. Here we observe that when 
d = 3 -|- 1 this effect is even stronger (due to the contribution of (1 -|- i^)/2) and in order 
to obtain a low temperature insulator it is sufficient to have 7 < \J{1 -\- i/)-^/[3(z^^ — 1)]. It 
is interesting to note that, when the normal phase is approximately insulating, the super- 
fluid phase transitions studied in section 4 are insulator /superconductor transitions. These 
are typically observed in the low temperature phase diagram of cuprate high-T supercon- 
ductors. At the same time we observe that when 7 = y^2/ (i/^ — 1) (for d = 2 -|- 1) and 
7 = (3 -|- v^l sj'iiy'^ — 1) (for d = 3 -|- 1) we have a linear in temperature resistivity, a 
phenomenon which is found in many unconventional superconductors^. 

^See also [46] for another mechanism to obtain a linear in temperature resistivity in holography. 
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In figure 3 we give the conductivity as a function of uj/T. That figure shows that 
Re[(T(ci;)] is large (small) when Z^^cp) is large (small) at the horizon; this behavior was 
expected because we showed above that it occurs at least for a; — )■ 0. 

5.2 Conductivity in the broken phase 

The broken phase should display an infinite DC conductivity as the system is supercon- 
ducting. We provide o" as a function of oj/T in figure 4. As w — )• the imaginary part 
diverges like l/u. This corresponds, through the Kramers-Kronig relation 

= -ip r <i.'S!i^ , ,5,15) 

■71" J-oo UJ'-UJ 

to a delta function in the real part: Re[o"(a;)] ~ irnsS^uj). This indeed means that the 
system is superconducting (a superfluid in the fluid mechanical case). The quantity can 
be identified with the superfluid density and from figure 4 we see that is generically 
changed by changing the couplings of the dilaton. 

Moreover, we see that the condensate leads to a gap ojg for charged excitations in the 
plot of Re[cr] versus cj, like in the S-BH configuration [11, 12]. Following [12], we identify 
Ug with the frequency at which Im[(T] has a bump (in figure 4 a minimum). We see that in 
the dilaton-BH, the gap depends on the couplings of the dilaton. We have also computed 
the ratio ujg/Tc to compare it with the BCS value, around 3.5, and with the result for the 
S-BH, around 8 with deviations of less than 8%, as discussed in [12]. For the values of the 
parameters chosen in figure 4 we find 

^ ~ 7.4 (8.7) for d = 2-hl and 7 = 0.4(-l) , (5.16) 

^-7.9(9.1) for d = 3-hl and 7 = 1(-1) ; (5.17) 

J-c 

therefore changing the couplings of the dilaton changes, although perhaps weakly, this ratio 
as well. The possibility of having such a high oOg/T^ compared to the BCS 3.5 is a nice 
prediction of holographic superconductor models: high-T superconducting materials show 
indeed a value which is systematically bigger than 3.5. We observe that a negative and large 
value of 7 corresponds to a bigger value of ujg (compared to that obtained with a positive 
7). It is interesting to observe that a large ojg is associated with an insulating normal phase 
according to the results of section 5.1. 

6 Inhomogeneous phases with broken U(l) symmetry 

All we have discussed so far can be applied both to superconductors and superfluids be- 
cause the dynamics of the gauge field was negligible: more precisely, since we are considering 
physical phenomena described by and O only, the backreaction of O on was not im- 
portant^. In this section, instead, we study inhomogeneous phases where the difference 

^However, from this statement we cannot conclude anything about the interaction between a^j and the 
condensate constituents, whose operators we have not introduced. 
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Figure 4. The conductivity as a function of the frequency in the broken phase, ^ > Hc, (solid bhie 
lines) and in the unbroken phase (dashed purple line). Re[(T] has also a term oc 6{uj) in the broken 
phase. We set T ~ 0.010/L, v = 1.2, Z^(0) = e'^'* and Z^, = 1. Moreover, we choose the bare 
kinetic term such that g'^ /el = L ln(z/L)|^=o for d = 3 + 1, while ei,/g — > oo for d = 2 + 1. 
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Generic quantity 


Superfluid (SF) 


Superconductor (SC) 




SF current density 


EM current density 


arg(C)) 


SF velocity potential in the lab frame 


condensate's phase 


ai 


external velocity in the lab frame 


EM vector potential 


in vortex phases 


SF 


SC 


aip-behavior 


a^p is frozen 


larger _ 

a<^ ~ n + ai^Jre ' 


tpci-behavior 


^cl ^ -000 (l - n^^VrM at B = 




quantization of B 


No 


yes: J dr rB = n 


free energy 


Fn - Fo n2 ln(i2/eGL) - nBR^/2 


finite as — )• c« 


first critical field 


largeR \ / r,2 n 


^ci = eg(i^i-Fo)/27r 


second critical field 


B,2 = 1/(2^2 J 


/7,2 = 1/(2^2 J 



Table 1. Comparison between superconductors and superfluids. A' and ^' are the penetration 
depth and coherence length; generically A' 7^ A and ^ ^, where A and ^ are the inverse masses 
of ai and ip^i = 1(0)1 respectively. ^, ai, i/'oo ^^nd are other constants. is the free energy 
per unit of volume Vd-s, of a vortex with winding number n and ^gl is the coherence length of the 
Ginzburg-Landau theory. The external magnetic field H is normalized such that it coincides with 
the total magnetic field B at 5* = 0. For superfluids we have H = B. See ref. [18] for a derivation 
of these model-independent properties using effective field theory methods. 

between a global and a local U(l) becomes crucial. In table 1 we provide a dictionary be- 
tween superconductor observables and the corresponding quantities in superfluidity, which 
are relevant for the following discussion. Given their importance in superconductivity and 
superfluidity, we study in particular vortex configurations. 

6.1 Superfluid vortex phase 

Let us start with superfluids. Even in this subsection we will stick to the superconductor 
notation in order not to introduce additional symbols; the reader can however translate to 
the terminology commonly used in superfluidity by means of table 1. 

The simplest example of vortex is a straight line, corresponding to the ansatz [21, 47, 48] 

* = V(^,r)e*"'^ , Ao = A^{z,r), A^ = A^{z,r), (6.1) 

where n is an integer. The straight vortex line is orthogonal to the plane spanned by the 
polar coordinates dy'^ = dr"^ + r'^dip'^, with < r < R and < < 27r, and passes through 
r = 0. We imposed r < R because, as shown in table 1, the superfluid vortex energy Fn — Fo 
diverges as i? — )• 00. We take R to be large compared to the radius size of the vortex core 
{i in table 1). 

Being irrotational, a superfluid does not rotate if the external angular velocity B/2, 
which we take to be constant for simplicity, is small enough. By using table 1 we see that 
B is related to the gauge field through 

= ^Br^ . (6.2) 
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The simplification which B = constant leads to is that it allows us to have a finite (J^) 
even for d > 2 + 1: we do not need to renormalize the divergence in (2.7) as d^J-yfj_ = 
in that case. So in this subsection we will also take e^/g — )• oo because we do not need 
any boundary kinetic term. For superconductors this will not be possible for d > 2 + 1 
as the dynamics of the gauge field forces B to be nonconstant (and in particular to decay 
exponentially for large r, table 1) as it will be discussed in the next subsection. 

Increasing B a single vortex appears and if B gets even bigger more and more vortices 
are created until a triangular lattice [49] is reached at S ~ Bc2, where i?c2/2 is the maximal 
angular velocity above which superfluidity is destroyed. The phase at i? ~ Bc2 is model- 
independent as any model of superfluidity reduces to the Ginzburg-Landau theory when the 
condensate becomes extremely small (which is the case at those high B). We therefore study 
the single vortex configuration here as it is energetically favorable when the holographic 
model is expected to deviate substantially from the Ginzburg-Landau predictions. 

Ansatz (6.1) leads to the following field equations 



W(^-^)/^fZ^{cl)) 



r 

+ 7* 



rf 



42 

p 



{A^ - nf 



= 



(A^ - n) = , (6.3) 



+ 



dr {rdr-Ao) 2WZ^{(P)ij'^ 



An 







{W^^--^y^ZA{(P)d,Ao) 
Wid-mZA{ct>) ' rf L^fZAi^) 

which have to be solved with appropriate boundary conditions. Our (external) conditions 
are (4.2) and (6.2) at the AdS-boundary, while at large r, r = R, we impose 

B 



(9r^ = , drAo = , A„ 



(6.4) 



The first two conditions correspond to the physical requirement that the field configura- 
tion goes to the homogeneous superfluid state, while the last one is a simple possibility 
compatible with eq. (6.2). Regularity requires the following boundary conditions at z = Zh 



f{zh)d,i: + 

f{zh)d,A^ + rdr (-drA 



drirdri)) {A^-nf 



4 = Q 



2WZ^{^)i;^' 
L^ZA[(t>) 



(A 



n] 



An 



(6.5) 



and at r = 



drAn = , A^ = Q 



— u , 



dri^ = for n = , 
'0 = for n 7^ . 



(6.6) 



We have solved eqs. (6.3) with the boundary conditions in (4.2), (6.2), (6.4), (6.5) and 
(6.6). In figure 5 we present the corresponding vortex profiles. 
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As shown in table 1, the maximum value of B for which the system is in the vortex 
phase, coincides with that for superconductors, that is Bc2 = Hc2- We have computed Hc2 
as the value of B at which ^' — )• and provide it as a function of /iL in figure 8 of the next 
subsection (dedicated to superconductors). 

6.2 Superconductor vortex phase 

In the superconductor case the dynamics of is important and to properly take into 
account we use the method described in section 2.1. In particular for vortices, the Neumann- 
like condition at the AdS-boundary, eq. (2.5), reduces to 



L 



d-3 



4 \r 



, (6.7) 



where we have set = 0. As already mentioned, the second term in the equation above 
is important in order to absorb UV divergences which occur for d > 2 + 1. Therefore, 
since we want to maintain the space-time dimension general, we have to keep that term. 
Moreover, the physical requirement that the field configuration goes to the homogeneous 
superconductor state tells us that at r — )• oo 

= , drAo = , A^ = n . (6.8) 

We have solved eqs. (6.3) with the boundary conditions in (4.2), (6.7), (6.8), (6.5) and 
(6.6). For n = the only solution we find is the homogeneous configuration of section 4 
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Figure 6. The condensate (left plots) and the magnetic field (right plots) for the n = 1 supercon- 
ductor vortex. Here we set v = 1.2, T ~ 0.15/L, ^^(0) = e"^*, Z.^,{(t)) = e""^ and et/g -> oo for 
d = 2 + 1 while we choose such that Cg ^ ~ 3.2L/g-^ for d = 3 + 1. 

with B = 0: this is the celebrated Meissner effect of superconductors. In figure 6 we give 
(O) and B = drUip/r as a function of r for the n = 1 vortex. Our solutions respect the 
model-independent properties given in table 1. This allows us in particular to obtain the 
coherence length ^' and penetration depth A'. We find 

C'c^ — , for d = 2 + 1(3 + 1) (6.9) 

for the values of the parameters chosen in figure 6, while we plot A' as a function of fiL 
in figure 7. Taking .^^(0) = e'^'^ for the sake of definiteness, we observe that when 7 is 
negative enough A' behaves in a qualitatively different way for large /i with respect to the 
S-BH. We observe that this is true both for d = 2 + 1, as already pointed out in [15], and 
d = 3 + 1. Since increasing /^iL for fixed TL corresponds to exploring the small T/Tc, region 
we conclude that the penetration depth in the dilaton-BH have a qualitatively different low 
temperature behavior. In particular A' does not increase more and more as T/Tc decreases, 
unlike in the S-BH [18]. We always find that A' diverges at fi = fic- this is a model- 
independent feature of superconductivity because when fi the condensate becomes 
small and the Ginzburg-Landau theory holds and predicts that divergence. 

Moreover, we have studied the first critical field Hd, the value of the external magnetic 
field H at the boundary between the homogeneous superfluid phase and the vortex phase. 
We can compute this quantity with the text book formula in table 1. We provide Hd 
vs. fiL in figure 8. Like for the penetration depth, in the high ^//ic (low T/Tc) region 
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Figure 7. The penetration depth as a function of fiL. Here we set i/ = 1.2, T ~ 0.010/L, 
Za{4>) ~ e^'^ , — 1 and Cb/g — > oo for d = 2 + 1 while we choose Cb such that ~ 0.41L/(7^ 

for c? = 3 + 1. We choose 7 = —4.4 on the left and 7 = 0.4 on the right. 

we find substantial qualitative differences compared to the S-BH case if the exponent 7 is 
below a certain bound: while for the S-BH (and for 7 above this bound) Hd — )• [18], 
in the other cases Hd remains sizable even at low T/Tc. This occurs because Hd oc Cq, 
whose low temperature behavior depends strongly on whether scale invariance is preserved, 
(f) = 0, or broken, (j) ^ 0- Working out eo explicitly one finds that in the scale invariant case 
Co — ?• and so does Hd, but when scale invariance is broken eo remains finite if 7 is below 
a certain bound. From eqs. (2.11) and (2.12) we see that this bound is —^^jiii^ — 1) for 
d = 2 + 1 and ^(v - l)2/(3(j^2 _ 1)) for d = 3 + 1. Interestingly, we observe that in both 
dimensionalities these values of 7 lead to an insulating normal phase at low temperatures, 
but the bound for d = 3 -|- 1 is substantially weaker than that for d= 2 + 1. This conclusion 
can be generalized to arbitrary ZpX^^ (not necessarily exponential): what actually matters 
here is whether ZA_^>)\z=Zh small enough for small T; but comparing with eq. (5.14), 
we see that when this occurs the normal phase tends to be insulating. Finally, the fact 
that Hd decreases for low [iL and eventually vanishes at /i = /ic is a model-independent 
feature of superconductivity: at the critical point both the superconductive and the vortex 
configurations go to the normal one. 

In the same figure 8 we also plot Hc2, the critical value of H above which superconduc- 
tivity is completely destroyed. We find that Hd < Hc2 for both d = 2 + 1 and d = 3 + 1. 
Therefore, the holographic superconductors under study here are of Type II, like their ana- 
logues without the dilaton [18, 27] and, interestingly, like all known high-T superconductors. 
The mere fact that we have Type II superconductors also implies that the phase at high 
external magnetic fields, H ~ Hc2, is a triangular lattice of vortices [49]. 

7 Outlook 

In this paper we have investigated various transitions in dilaton holography, including in 
particular those that break a local (superconductor) or global (superfluid) U(l) symmetry 
in diverse dimensions. A detailed summary of the results found is provided in the second 
part of section 1. Therefore, we focus here on possible outlook of the present work. 

In this article we have taken the limit in which Gn ^ so that the matter fields (the 
gauge field and the charged scalar) do not change the metric and the dilaton. The next 
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Figure 8. The critical magnetic fields as a function of (the inset is an enlarged version of 
the plot of Hci) . Here we set v = 1.2, T ~ 0.010/i, Za[(\)) = e'^'^, 7 = -4.4, Z^{(j)) = 1 and 
eh/g — > oo for d = 2 + 1 while we choose Cf, such that e^'^ ~ OAlL/g^ for d = 3 + 1. 



interesting step would be to include the backreaction, that is going beyond the Gat — )• 
limit. A first step towards this goal could be analyzing the homogeneous superfluid phase of 
section 4 away from such limit. Then, a second (more ambitious) goal might be finding the 
backreaction of the vortex solutions on the dilaton-gravity system. Since a vortex carries 
angular momentum, one would expect the resulting bulk solution to resemble a rotating 
charged black hole. 

Another interesting extension of this work is the quest for a unified holographic model 
for the full cuprate phase diagram which includes the insulating (ferromagnetic), supercon- 
ducting, strange metallic [50] and pseudogap regions. 
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